Specific Heat of Quantum Elastic Systems Pinned by Disorder 



^1- 
O 

o 

> 
O 

OS 



c 

o 
o 



> 

(N 



^1- 
o 

-I— > 



c 

o 
o 



X 



Gregory Schehr, 1 Thierry Giamarchi, 2 and Pierre Le Doussal 3 

1 Theoretische Physik Universitat des Saarlandes 66041 Saarbriicken Germany 
2 DPMC, University of Geneva, 24 Quai Ernest- Ansermet, CH-1211 Geneva, Switzerland 
3 CNRS-Laboratoire de Physique Theorique de I'Ecole Normale Superieure, 24 Rue Lhomond 75231 Paris, France 

(Dated: February 2, 2008) 

We present the detailed study of the thermodynamics of vibrational modes in disordered elastic 
systems such as the Bragg glass phase of lattices pinned by quenched impurities. Our study and 
our results are valid within the (mean field) replica Gaussian variational method. We obtain an 
expression for the internal energy in the quantum regime as a function of the saddle point solution, 
which is then expanded in powers of ft at low temperature T. In the calculation of the specific heat 
C v a non trivial cancellation of the term linear in T occurs, explicitly checked to second order in 
ft. The final result is C v oc T 3 at low temperatures in dimension three and two. The prefactor is 
controlled by the pinning length. This result is discussed in connection with other analytical or 
numerical studies. 

PACS numbers: 



I. INTRODUCTION 

Despite progress in analytical solutions of models of 
glasses in solvable limits such as mean field, there is at 
present little known information about detailed thermo- 
dynamic properties of experimental interest such as the 
specific heat. In particular an outstanding question is 
its linear behavior at low temperature 1 ". Such behavior 
was measured in a large variety of experimental systems 
including amorphous solids, disordered crystals and spin 
glasses 2,3,4 ;. A phenomenological approach, based on the 
existence of two level systems^ was proposed leading to a 
linear temperature dependence of the specific heat. De- 
spite the success of such prediction for many experimen- 
tal systems, the question of the validity and applicability 
of such arguments is still under debate. In particular, 
for many glassy systems, the microscopic origin of the 
assumed two level systems remains unclear. 

Computing the specific heat from a microscopic model 
is of course an extremely difficult problem for a dis- 
ordered system. However, mean field theory of quan- 
tum spin glass, have been studied^ leading to results 
in agreement with the linear dependence of the specific 
heat. Nevertheless the validity of the linear tempera- 
ture dependence is being challenged^. Another impor- 
tant class of glasses to which these mean field methods 
have been applied with success to compute correlation 
functions 1 ^, consists in disordered elastic systems. Such 
systems cover a wide range of experimental situations 
such as charge density waves^i, electron glasses 1 ^ 3 -, and 
flux latticesi 4 * 1 ^* 1 ^*^, for which the quantum limit is of 
interest. Thus both from the experimental and theoreti- 
cal point of view such disordered elastic systems are ideal 
to address this important question of the behavior of the 
specific heat. 

This is the question that we address in the present pa- 
per. We consider these systems in the elastic limit when- 
ever topological defects can be neglected. This is the 
case within the Bragg glass phase which was studied pre- 



viously in both classical^ and quantum 1 ^ limit using the 
Gaussian variational approximationi2ii2iiS to the repli- 
cated Hamiltonian. We use in this paper the same vari- 
ational approach to compute the specific heat for these 
systems. In two previous short papersSiSi, we have sum- 
marized the result of the calculation of the specific heat 
for these systems in classical and quantum limits, as well 
as the applications to superconductors. In the present 
paper we give a detailed exposition of the method. In- 
deed two important questions must be addressed. In or- 
der to obtain the correct temperature dependence, it is 
necessary to push the understanding of the structure of 
the mean field solution beyond what has been achieved 
in Ref^. In addition obtaining the amplitude requires 
performing a semiclassical expansion. Besides, for sys- 
tems in dimension d < 2, two saddle points solutions 
for the variational equations exist, one coming from the 
thermodynamics and one from the so-called marginality 
condition. We discuss the respective roles of these solu- 
tions. 



The outline of the paper is as follows. In Sec. ITT1 we 
formulate the problem of the calculation of the internal 
energy for an elastic system with quenched disorder. In 
Sec. IIIII we recall the replica variational method and ob- 
tain a compact expression for the internal energy in terms 
of the saddle point solution. This expression being very 
hard to compute analytically we present a systematic ex- 
pansion in powers of h, given that the specific heat is a 
function C v (f3h,h). In Sec. IIVI we compute the first two 
terms of this expansion, the details being given in Ap- 
pendix ^ and Appendix [BJ The result is a cancellation 
of the linear term, resulting in C v ~ T 3 . 
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II. SPECIFIC HEAT COMPUTATION USING 
FUNCTIONAL INTEGRAL. 

A. Model 

We consider a collection of interacting quantum ob- 
jects of mass M whose position variables are denoted 
u a {Ri,r). The equilibrium positions Ri in the absence 
of any fluctuations form a perfect lattice of spacing a. 
Interactions result in an elastic tensor <I> Qj ^(g) which 
describes the energy associated to small displacements, 
the phonon degrees of freedom. Impurity disorder is 
modelled by a t independent gaussian random poten- 
tial U(x) directly coupled to the local density p(x) — 
J2i $( x — Ri~ u(Ri,r)). We will describe systems in the 
weak disorder regime a/R a <C 1 where R a is the trans- 
lational correlation length, e.g. in a Bragg glass phase 
where the condition \u a {Ri, t) — u a (Ri + a, r)| <C a holds, 
no dislocation being present. The system at equilibrium 
at temperature T — 1/(3 is described by the partition 
function Z = Tre- [iH ^ n ^ = j DuDU e - s ^ h with the 
Hamiltonian H [II, u] = i? p h[n, u] + i?di S [n, u]: 

1 f H(a) 2 

ff P h[n,u] = - J -^ + ^« a (# ai( j( 9 m- ? ) 



a, 



ffdisP, «] 



d d xU(x)p(x 1 u(x)) 



(1) 



/ dr in Q {q,T)d T u a {q,T) + H[Il,u} (2) 

JQ J a 



and its associated Euclidean quantum action in imagi- 
nary time r 

r/3fi 

5fn,ui = - 

/0 Jq 

where u(q, r) and its conjugated momentum H(q, r) sat- 
isfy periodic boundary conditions, of period j3h, along 
the r axis. One denotes by J q = J BZ j§^jd integration on 
the first Brillouin zone. For simplicity we illustrate the 
calculation on a isotropic system with $> a .f3(q) = cq 2 8 a j3 



and denote disorder correlations U(x)U(x') = A(x- 
with U(x) = 0. 



■ x 



with f2 the volume of the system. _ff p h [II, u] being 
independent of the imaginary time r, Hp^Il, u](t) — 
e TH pt [n,, l ] ffphe -rH ph [n, u ] = H ph [U,u], one writes 
(ffphp, u]) as 



(H ph [U,u}) =Z pX ] / VXWue 



1 



drH ph [II, u](t) 

where S p h is the Euclidean quantum action © corre- 
sponding to 10. One introduces the Fourier transforms 
(w.r.t. imaginary time variable) of the fields in terms 
of Matsubara frequencies to n = 2-nn / (0h) , u(q,u> n ) = 
Jy h dre tUnT u(q, r) and similarly for II(q, u> n ). After inte- 
gration over the field II in J5J , one obtains by performing 
the change of variable H(q, r) — * H(q, r) + iMd T u(q, r)): 



S <flphp,tt]) 



Of 



Musi 



cq z 



(6) 



Performing the sum over Matsubara frequencies, e.g. by 
using the spectral representation of the 2-point Green 
function l/(Mw„ + cq 2 ) one obtains 



prZP 



Kv / \q\f B (v\q\) 



(7) 



where is the (temperature independent) zero-point 
energy, /b(x) — (e^ hx — 1) _1 is the Bose factor and v = 
y/c/M the pure phonon velocity. Computing the specific 
heat (0} from Q) one recovers the Debye law for pure 
phonons 



C v Bebye(T) — Ad 



I_\ +0 (T d+l ) (8) 
nv ' 



where A& = T(2 + d)C(l + d)Kd, with £(x) the Riemann 
zeta function, Kd = Sd/(2n) d , Sd being the volume of 
the d-dimensional unit sphere, e.g. A 3 — 2tt 2 /15. This 
method using a functional integral formulation (jSJ) is thus 
a very convenient way to compute the specific heat of 
elastic systems. 



B. Pure phonons 



C. Disordered case : general expression 



We first consider the case of pure phonons, described 
by a purely isotropic elastic hamiltonian, i.e. with 
U = 0: 

ffphP, U] = I [ + cq 2 u(q)u*{q) (3) 



M 



One computes the specific heat per unit volume C v (T) of 
this system using the functional integral in imaginary 
time ©, C V {T) being defined by 



C V (T) 



1 d(H ph [U,u] 



n dT 

(H ph [Tl,u]) = -LTr[i/ ph [n, U ]e-^[ n '"]" 



(4) 



We now extend this analysis to the disordered case Q 
and obtain the analogous formula of ijfJJl for the disordered 
average internal energy. We remind the expression of the 
disordered elastic hamiltonian 

#[11, u] = H ph [U, u] + [ U(x)p(u(x), x) (9) 

U(x)U(x') = A(x-x') 

and denote by P(U{x)) the gaussian distribution of the 
disorder. As p{u{x),x) depends only on the field u(x) and 
the disorder is r — independent, we compute (H[H,u]) 
in one realization of the disordered potential in the 
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same way as for the pure case (0 with the substitution 
iJpiJII, u] — > H[H, u] ©■ To compute the average over 



the disorder, we introduce replicas to deal with the de- 
nominator in (JSJ, using Z^ 1 = lim^o Z • 



(-ff[n,d) = lim f VUPiUix)) [ VIi a Vu a ^- f drH[Il ai , u ai ]{T)e-^ £<» s[n 



(10) 



where 5 [11°, u a ] is the action defined by and a — 1, ..k 
is a replica index, a± being one of these replicas, and 
VU a = PIT 1 ... IF and similarly for Vu a . As for the pure 
case, the gaussian integrals over the fields IP is straight- 
forwardly computed (performing the same change of vari- 
able). To treat the integral over the disorder U(x), one 
uses the identity for a gaussian variables with correlations 
A(x-x') © 



/ : C7exp(-/ : U) = -f : A : /exp ( -/ : A : / ) (11) 



for any vector / = f(x) and where / : U = 
.(, f(x)U(x) and / : A : / = f , f(x)A(x - 
x')f (x'). Applying this formula (fTTfl to Qjjf with f(x) = 

I J2a fo^ drp(u a (x, t), x) ©, one obtains the exact for- 
mula for the internal energy: 



A(x - x')(p(u a (x, T),x)p(u b (x', r'),x% ep 



dr j dr' 

^0 -ix i j 



,6=1 



(12) 



where G(q,uj n ) = (u a (q, uj n )u a *(q, w n )) Yep and the aver- The same kind of manipulations 10 lead to the following 



aged values in ((TUl (...) rep = / Vu a ...e ~ are com- 
puted with the replicated action given below l|15fl . Writ- 
ing the density p(u(x),x) asiS 



p(u(x), x) = p (l - V • u{x) + e iK < x - u ^) (13) 



replicated action 



ST*[u] = fj2 i( c ? 2 + Mw 2 n )u a (q,UJ n )u a (-q,-UJ n ) 
J 1 n,a Lp(l 

ddxdTdT ' ZZ R ( U "( X > T ) - U b( X > T ')) ( 15 ) 



In these expressions iflU fT5)l . Ak = J d d xe lK ' x A{x) de- 
note the harmonics of the disorder correlator at the re- 
ciprocal lattice vectors K . The exact expression l|14|l is 
with K the reciprocal lattice vectors of the perfect lattice the star ting point of our computation of the specific heat, 
and po oc a~ 2 the average density, the second term of fl^jl 
can be written as, discarding irrelevant terms: 

III. VARIATIONAL COMPUTATION 



A- 



Given the complexity of the replicated action 
l)15|l. we study it within the Gaussian variational 
approximation's^. The use of such a method for the cal- 
1 l r l JL, / culation of the specific heat necessitates to compute the 

fc/3ft2 J d X ^ T ^ T (R( u a( x i T ) — u b(x,r )))rep full temperature and ft dependence of the saddle point 

solution. Before establishing the formulas for the specific 
(14) heat in Sec.llTTGl we present the full variational solution 
in Sec. ESI and Sec. UTTE! 



a, 6=1 



R ( u ) = Po zZ Ak C0S ( K ' 

K 
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A. Saddle point equations 

The variational methodifl*^ is implemented by choos- 
ing a gaussian variational action So, parameterized by a 
k x k matrix in replica space G~£(q,u) n ): 

S ° = Vrp, f ^2 G ab( ( l> UJ n)u a (q,ui n )u b (-q,-ui n ) 

L P nJ 1 a ,b 

G~l{q,uj n ) = cq 2 S ab - a ab (16) 
which minimizes the variational free energy 



F™ = F + —(S^-So)so 
Fa = \j J2(^G) aa (q,u n ) 



(17) 



The disordered term in S lcp 1)15)) being purely local in 
space, but bi-local in time the self energy does not depend 
on q and depends only on oj n , thus a a b = CatK)- How- 
ever, the disorder potential being r-independent, non di- 
agonal terms (in replica space) such that o a =£b do exist 
only 18 for uj n = . Indeed, before averaging over the 
disorder, different replicas are independent: 

GW = (u a (x,T)u b (0,Q)) = (u a (x,T))(u b (0,Q)) 
= {u a (x,0)){u\0,0)) (18) 

In the limit k — > 0, one parameterizes G aa {q, w n ) by 
G(q,tu n ) and G a =£b(q,w n ) par G(q,u), < U < 1, and 
in a similar way B a ^{r) — ([u a (x,r) — u b (x,0)] 2 )) /N 
by B(t) and B(u) independently of r. Using the inver- 
sion rules of hierarchical matrices one obtains the saddle 
point equations (directly written in the limit k — > 0) by 
minimizing the variational free energy F™ (|17l) . 

G-^q^n) =J2 G ab(Q^n) = cq 2 + Mif n (19) 
b 

2 fP h f 1 
+ - dr(l-cos(w„T))(F'(S(r))- / duV'(B(u))) 
11 Jo Jo 

a{u,uj n ) = Snfi^^V' '{B(u)) 



with 



2h f 



+ S I ^ (<Z ' W " = 0) ~ G{<1 ' Un = °' U) ) (20) 
os r 

^ = Aft E / - cos ^« t )) ( 21 ) 

with V(B) = -/OqEk A K exp(- J Bi4: 2 /2). An interest- 
ing property of the above saddle point l|19(l is thus that 
replica symmetry breaking is confined^ to the mode 
u) n = 0. The equation for this mode 



a(u,u) n ) = 5 nt o?j^V'(B{u)) 



has been studied previouslyi&( notice that this is identi- 
cal for the one for a model with point-like disorder in 
d dimensions studied inifi^). For the potentials with 
power law correlators V(x) = gx 1-1 / '(2(1 — 7)) there are 
two generic cases: long range correlations 7(1 — 2/rf) < 1 
for which one has full replica symmetry breaking (RSB) 
and short range correlations 7(1 — 2/rf) > 1, for which 
one has a one step RSB (and a transition to high tem- 
perature replica symmetric (RS) phase). For the single 
cosine model defined by 1)14(1 where the sum over K in 
R(u) is restricted to the lowest harmonic, there is a one 
step RSB solution for d < 2 and a full RSB solution for 
d > 2. More realistic models involving several length 
scales include several of the above regimes^. The previ- 
ous analysis of this equation 1)22)1 revealed the existence 
of a breakpoint u c such that cr(u) — cr(u c ) for u > u c . 
In the case of a full RSB solution, cr(u) is a continuously 
varying function of u for u < u c . In d < 2, the single co- 
sine model is instead described, in the low temperature 
phase, by a one step RSB solution such that cr(u) = for 
u < u c . 

Independently of the RSB scheme, one can write the 
variational equations 

G~ 1 (q,u n ) = cq 2 + Miu 2 n + E(l - 6 n ,o) + I(u n ) (23) 
2 ,/J» 



I(u n ) = lJ o dr(l - cos (w„t))(F'(S(t)) - V(B)) 



B = B(u>u c ) = -Y J j 



1 



cq 2 + Mul + E + JK) 



together with the expression of -B(r) 121)) and where we 
have used the definitions E = [<t](u c ), [c](u) = ua(u) — 
J dva(v) together with the equation 1)22)1 . 



B. Solution of the variational equations 

The general way to study the equation (|2^)l has been 
presented in details ini&. Here we present only the main 
results relevant for our study. 



1. Periodic structures in d > 2 : marginally stable solution. 

In that case there exists a full RSB solution, the func- 
tions B(u), a(u) and [<r](u) being obtained by elimination 
from the system 



1 = -4V"(B(u))J 2 {[a](u)) 
<u) = 2 -fv'(B(u)) 



(24) 



Jn{x) 



1 



q {cq 2 + x) r 



(22) 



Once this system ()24)l is solved for u < u c , the constants 
E and B are unambiguously determined by Eq. 1(23)) 



5 



together with the so-called marginality condition 



1 



-4V"(B)J 2 (E) 



S (4-d)/2 = ^LV"{B) 



Old = 



(25) 
T vK d (d - 2) 
sin 7rrf/2 



where the second line in (|25|l is valid only in the infinite 
UV cut-off limit. In order to understand the finite tem- 
perature behavior of the variational equations below u c , 
it is useful to write the equations (|24|l in terms of the 
rescaled variable w = BHu/h. One thus has 



Bhu 
W = — 

Bh 

<7{U) = —S(W) 

[a](u) = [s](w) 



(26) 

(27) 
(28) 



where the function [s](u>) = ws(w) — dvs(v) is h and 
(3h independent. Indeed, differentiating (|24|l . one obtains 
that it is implicitly defined through: 

w = 4 ( ^ ([ " ])) V w ((-^)- 1 ( t ^ 1 t)) (29) 

Notice that for the single cosine model corre- 
sponding to V(B) = — Wexp(— K 2 B/2) one has 
V»'((-V»)--L( x )) = K 2 x/2, and Eq. @ gives back 
[s]( w ) = ((4 - d)w/{2K 2 c d c- d ' 2 )) 2 ^ d -^ with c d = 
(2-d)^ 1 -^/(2 d + 1 sin(7rd/2)r[d/2]) {e.g. c 3 = l/(8n))HL. 
The function [er](u) = [s](w) being independent of H and 
BH iP?|). it follows from that B ( u ) = B i w ) is also 
independent of h and (3h below the breakpoint u < u c . 
The equation l|29|l written at the breakpoint gives 



u> c (E) 



Bhu c 



= * { -^§fv>"(B) (30) 



Finally, once the equation l|29(l is solved, w c , S, I(ui n ) 
which all depend both on h and H are determined by 
Eq. lO, H3) and ^ together with the definition (ED- 
Notice also that combining (|25|l and 130(1 . one obtains 



w c 6B + 2J 2 (E)5Z = 



(31) 



where S stands for an infinitesimal variation : a useful 
identity, only valid for a full RSB solution, in the follow- 
ing computations. 



2. Periodic structures in d < 2 
margino 



thermodynamics vs 



The case d < 2 is more subtle, since as noticed 
previously^ the saddle point equations admit two so- 
lutions. Which solution to choose was shown to be im- 
portant for the transport properties and will of course 
be important for the specific heat as well as will be dis- 
cussed in Sec. lIVBl Let us examine here the two possible 
solutions. 



For the single cosine model corresponding to V{B) = 
— Wexp(—K 2 B/2), in dimension d < 2, the variational 
equation (|2*2"|l admits a one step RSB solution given by 10 



[a] (it) = , u <u c 



[a](u) 



Bhu c -V'{B) 



(32) 

u>u c (33) 



In this one step RSB case, one can use only (|32l I23|) 
to determine the three quantities B, £ and u c : thus 
one equation is missing. In the statics, the breakpoint is 
then usually obtained by minimizing the variational free 
energy F val ' with respect to u c , the so called thermody- 
namical condition. It leads to in d < 



Bhu. 



th 



^ w 



th 



K 2 h^M^ h ) 



(34) 



where the second line in l|34|l is valid in the infinite UV 
cut-off limit with a d given in l|25|l . However, it is known 
that this condition (|34|) gives incorrectly the behavior 
of dynamical quantities such that the conductivity^. A 
distinct choice is to impose the marginality condition 125|) 
which, given the relation l|32l) 



-4V r "( J B mg )J 2 (E mg ) 

= _ K 2^d J: (d-2)/2 c -d/2 

4 



(35) 



which allows to obtain the correct dynamical behavior. 
One can show, using a Keldysh mean field approach and 
performing analytical continuation to imaginary timaS^, 
that this is indeed the correct solution from the dynam- 
ical point of view, i.e., if one considers, in an infinite 
system, the large time limit where time translational in- 
variance and equilibrium fluctuation dissipation theorem 
hold. 

Finally, for the special case d = 2, these two conditions 
(|34|l and ll-iol) coincide and there is no ambiguity in that 
case. Therefore, in the following we will treat this case 
together with the full RSB one in d > 2. 



C. Internal energy (H) within variational method. 

In this Section we use the variational method, and 
identities valid at the saddle point, to derive a compact 
and useful expression for the internal energy (H) <|14|) 
which is analyzed in the following Sections. The idea 
is indeed to compute the averaged values in l|14f) with 
the trial gaussian action Sq (|I6fl instead of the exact one 
5 rcp l|15fl . One can show that thanks to the variational 
equations, it is equivalent to compute C V (T) using the 
variational free energy C V (T) = -Td 2 F val /dT 2 in- 
stead of the exact one. 

We start by deriving some identities which will be use- 
ful below to replace the kinetic term cx j q J2 n c 1 2 G{q, w„) 
in l|14|) by a more convenient one. First we remind that 
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G(q,u> n ) — G c (q,uj n ) for u> n 7^ (a consequence of the Combining 1(371 138(1 and using the definition of B(u) 1231) 
r-independence of the disorder dJ)). Using simply the one obtains 
definition of <r a b(q, ui n ) i|16J> an( i G~ 1 (q, iv n ) 1)19(1 one has 



<?aa (<7, 0J n ) = cq 2 - G c 1 (q, U) n ) - 5 n>0 ^ a ai> ( 36 ) 

Using (pTfTf) together with the identity 

SbG' ab (g,w„)G' fcc 1 (g, uj n ) = 5 ac , one obtains in the 
limit k — > 

cq 2 G{q,uj n ) = cq 2 G c {q,uj n ) (37) 

+ S n,0^Y J dU fa q > W « = °) - G (9- «)) 

where we have used G~ 1 (q,uj n = 0) = eg 2 and the vari- 
ational equation 1)22(1 . Using the variational equation for 
I(w n ) and the identity J* duV'B(u) = V'(B) - one 
has 



^2 / c 9 2 G(g,w„) 

n J 1 

/" C g 2 + S + /(W„) 



ft 

ft 



(39) 



2/3H ^J q cq 2 + Z + MojI + I(w„) 



2ft 



/3ft 



— / drfl(r)V'(B(T)) + £r / dui?( U )U'0B(u)) 



2// 



— £ G c ( 9 , Wn )(S(l - + J(w n )) = 



(38) 



This last expression ((3911 allows to write a compact ex- 
, pP% „i pression for (H) (114(1 computed within the variational 

— / (1tB(t)V'(B{t)) - V / dttG c (g,w„)V' / (B(u)) method under the form 

I 



eg 2 + S + J(gn) _ 1 

eg 2 + £ + Mw 2 + I(u n ) + ft 



0h 



cIt[F(B(t)) — F(B)] — / dw[F(B(w)) — F(B)] (40) 



r 



where F(X) = V(X) - f V'(X) and w c given in ffity. 
Although this form is compact and convenient, its tem- 
perature dependence is hard to extract, and we will resort 
to an expansion in powers of ft. 



heat to lowest order. 

In order to check whether the obtained temperature 
dependence is correctly captured by the lowest order cal- 
culation, we examine in Sec. the higher orders. 



IV. SEMI-CLASSICAL EXPANSION: LOWEST 
ORDER 

In this section we extract the temperature dependence 
of the specific heat from l(4U|) . In order to do so we use an 
interesting property of the variational equations ((23125(1 : 
the solution, as well as the internal energy 1(40(1 can be 
organized in an expansion in ft keeping /3ft fixed. To 
tackle analytically these variational equations, we will 
thus organize our calculations using this expansion for 
any quantity (not designated by a calligraphic letter) 
Q(ft,/3ft)=££°ft"Q n (/3ft). 

The calculation of the specific heat to a given order 
in the expansion requires the knowledge of the saddle 
point solution quantities, including I(ui n ), to the same 
order. Such an analysis of the saddle point solution was 
performed previously^ only to lowest order. We start by 
recalling this analysis and extracting from it the specific 



A. d > 2 

The analysis of the variational equations 1(231 125|) leads 
to the following equations for Io(u> n ) 



JoK) = -4U"(0) (ji(S ) - Ji(S + + IqM) 

(41) 

and for Eq 



1 = -4V"(0)J- 2 (£ ) 



(42) 



Because of the marginality condition 1(42(1. the function 
Io{w n ) is non analytic, its low frequency behavior being 
given by 



/oK) ~ J^flJyKI +°("D ( 43 ) 
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The marginality condition 1)42)1 thus leads to a gapless 
excitation spectrum. Indeed, the low frequency behav- 
ior of the analytical continuation Io(u) n —iu + + ) = 
Iq(ui) + il{f(u>) reads 



I {lo) ~Alu 2 



-Blj 



(44) 



with A=%{1- &jjf) and B = where J n = 

j7n(So)- Notice also that at this lowest order the equa- 
tions <f 4 11 142)1 show that Io(uj n ) and Eo are independent 
of [3H (Io(u) n )) depends of course implicitly of 0H through 
Lo n ). One then obtains the lowest order expansion from 

601 {H)/n = h q + hHxipn) + o(h 2 ) 



Hn 



-2F'(0)J X (S ) 



dw[F(B(w)) - F(0)] 



1 



eg 2 + E + I {uj n ) 

q Cq 2 + -^UJ 2 + E + loi^n) 



-F'(0)[2S 1 J 2 (E ) + U ; COJ B 1 ] 



(45) 



where we have used that B(w) is ^-independent for w < 
w c as well as lim^,^,^- B(w) = B in the case of a full 
RSB solution and F"(0) = 0. From JSTJ, the last terms 
in expression l|45|l just cancel leading simply to 



H = 



1 



eg 2 + E + I (u n ) 

cq 2 + AujI + E + JoK) 



(46) 



From Eq. (J20J, one has that w c0 = w c0 (So) is inde- 
pendent of (3h, as well as B(w) for w < w c such that 
Hq is (3h independent. To compute the specific heat to 
lowest order in this h expansion, C v = C v0 (/3h) + 0(ti), 
one thus focuses on Hi, whose temperature dependence 
is contained in the Matsubara frequencies. Transforming 
the discrete sum over Mastsubara frequencies in l)4(i)l in 
an integral one obtains, 



/ + O0 T 
-oa T 



p(ui) = -rw / lmG c (q, UJ n — > —ioJ + + ) 

V Jq 



-iM 



(47) 



£48) 



V 2 ( C g2-^ + So+J , M 2 + (W) 2) 



where p(uj) is the density of states. All the temperature 
dependence in 147(1 is now contained in the Bose factor. 
It follows from g7J| that 



C V0 (T) = (f3hf 



dto p(u>)lu 2 
4tt sinh 2 (3hw/2 



(49) 



The low temperature behavior of the integral in H49|) 
is governed by the low frequency behavior of the den- 



sity of states p{u>) 



u> 



(|4*51 ISJ, which then leads to 



Cvo(T) ~ T 3 in all dimensions d > 2. To this lowest 
order, there is thus no linear nor quadratic term in T in 



the specific heat. The specific heat has the dimension 
of the inverse of a volume, its low temperature behav- 
ior is then in general characterized by a typical length 
scale and an energy scale. In our problem, a natural en- 
ergy scale T p is given by the pinning frequency uj p with 
T p = fkUp = hv^Yjo/c — liv/R c , where R c is the Larkin 
length. From 1)490. one obtains 

C v (T)~^-K d R- d F Cv [Rc/a] ^ ' ™* 
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-0(h, [T/Kf) 
. ' ::,,M 

where Fc v (x) is a scaling function with asymptotic be- 
haviors given by 



FcM 



FcM 



1 



V4^d 
2 d+i 



d-2 



sin nd/2 



x>l (51) 
x < 1 (52) 



Thus, at very weak disorder, i? c 3> a 1)5 ip. and the typi- 
cal volume associated to the specific heat is R d , although 
at stronger disorder R c <C a (I52|l it becomes a d , the scal- 
ing function Fc v (x) describing the crossover from one 
regime to the other. Notice also that, in d = 3, in the 
limiting case R c /a — > oo, l|5U|) together with l|51|) and 
K3 = l/(2-7r 2 ) give back exactly the Debye formula © 
for pure phonons. 

Our results are thus at variance with the naive linear 
T expected from the two level system phenomenology. 
Contrarily to the pure case, the specific heat is now pro- 
portional to T 3 independently of the dimension. This 
specific heat reflects the density of states of the modes 
of vibration of the disordered system. A formula such 
as (|49f) shows that the specific heat would be the same 
than the one for free independent vibrational modes with 



a density of states p{u>) 



Note however, that since 



here the modes are not independent the precise calcu- 
lation of the specific heat cannot be done naively and 
requires the full calculation of the energy, as we have 
performed in the present paper. 



B. Periodic structure in d < 2. 

In that case, where there is a one-step RSB solution, 
the equation for J (w n ) is also given by (|4*T1) but the 
choice of the thermodynamical condition 1)34)1 leads toi£ 
/j h (w„) cx lu 2 , for small v n , and indicates a gap in the 
energy spectrum. 

Instead, if one uses the marginality condition i)35)l . the 
behavior is then similar to the previous full RSB solution 
within I™ 6 (uJn) oc \u) n \, an d describes a gapless energy 
spectrum. As discussed in the previous section IIII B 21 
the marginality condition is the consistent prescription 
(within this Mean Field approach) to compute dynamical 
quantities. 

We now turn to the computation of the internal energy 
in both cases 1)341 135)) . In that case, where there is a one 
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step RSB solution, the function B(w) is now discontinu- 
ous at the breakpoint w c . From this discontinuity results 
an other contribution to H\ compared to expression for 
the full RSB case (|45|l . Indeed in that case one has in 
both cases: 

H = -2F'(0)Ji(E ) + w cO F{0) 

H _ J_ f eg 2 + S + I (uj n ) 

+ J F'(0)[2S ll 7 2 (E ) + wM + w cl F{0) (53) 

where Hq is still temperature independent. Notice that 
here, one can not use the general formula l|31ll valid 
only for a full RSB. If one uses the thermodynami- 
cal condition, one has for the periodic case V(B) = 
-Wexp{-K 2 B/2) 



Hi 



m ^ J q cq 2 + ^loI + S* h + I*(u n ) 



AH! 



th 



AiJ! th = V(0) 



„th 



Si 

2S^ C 



th 



K' 2 



£i h j2(i^ h ) 



carefully to distinguish between thermodynamics and dy- 
namics quantities occurs also in other physical quantities. 
The capacitance (or compressibility) is also different de- 
pending on whether one considers the thermodynamics 
one or the small frequency onaiSii^iS^. The calculation of 
the specific heat clearly shows that for a thermodynamic 
quantity one should not use the criterion given by the 
marginality condition. The discussion of Sec. lIII B~2l sug- 
gests that the specific heat defined by the marginality 
condition (i.e. negative in the present case) could cor- 
respond to an experiment on a aging system using slow 
time dependent heat. 

Using the thermodynamic saddle point equation, on 
the other hand gives in one dimension (more generally for 
d < 2) an exponentially small specific heat at low tem- 
perature due to a gap. This is clearly an artefact of the 
variational approach. Corrections away from mean-field 
will most likely transform this gap into a pseudo gap, and 
yield a T a behavior for the specific heat. It is reasonable 
to surmise that the exponent a is larger than three and 
thus that the variational method attempts to reproduce 
such a large power by artificially inducing a gap. Indeed 



some studies of one-dimensional system: 



,24.25,26 



suggest 



4 "c u E th 



^0 



8 sin7r<i/2 



(54) 



where we have used <|34[) . One thus obtains formally the 
same expression as previously obtained for the full RSB 
case (@I)J). But here, as 7o h (w n ) cx w 2 , describing a gapped 
excitation spectrum, one obtains that the specific heat 
C* h (T) vanishes exponentially at low T. If one imposes 
instead the marginality condition Eq. I|35fl one obtains, 
using 



H™ s 



P h n Jq cq2 



c^ + s^ + c^K) 



a power law behavior p(w) ~ to . Such a density of state 
would lead to a T 5 behavior for the specific heat. 



V. BEYOND THE LEADING ORDER 

Given the fact that the linear temperature dependence 
disappears in (|49H . it is important to know whether this 
property holds to higher orders. We thus derive explicitly 
the solution to the next order. In doing so we unravel a 
structure of the variational solution not present to lead- 
ing order. 



A#i mg 



AiJ! mg = 



2 \2 V"{0) 2 
1 (2V(Q)V"'(Q) V'(0)V"'(0) 



4-dy V'(0)V"(0) V"(0) 2 

rpp. d 2 H v; "V 



(55) 



A. A need for a singular self energy in d > 2. 

We first need to extend the study of the variational 
equations to next order, allowing to extract Si and 
Ii(u> n ). To this order, as we will see, this quantity be- 
comes (3h dependent. From (|23l) and l|25l) . one obtains 



J ° 2(A-d)(3h^ 



r? c^+^+sr+/ ras K) 

This results in the low temperature behavior of C™ g (T) 



Si = 



V"'(0) 



-8(U"(0)) 2 J3(S ) 



C™ S (T) = (f3h) 2 



Sl + h(u> n ) 



1 + 4V"(0)/C 2 (a; n ) 



BlI (uJr, 



no) 

1 V"(0) 



sinh 2 j3Hu/2 



^ d-2 t 



(56) 



W"'{0) 

+ pn 



y]/Ci(a; m )(fci(a; m ) - £i(a; m + u n )) (57) 



(58) 



which is linear in T but negative for d < 2. 

This shows explicitly in that case the inconsistency 
of the marginality condition (|25|l to compute thermody- 
namical quantities. That for such glassy system one has 



, (cg 2 + M W 2 +S + /o(0)P 



Note that in the limit ui n — > the denominator in the 
l.h.s of (123 behaves as -8V" (0)/ K0^3( s o) and thus 
the first term yields exactly Si in this limit. Given the 
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complexity of the second term in the l.h.s of i|57|l . we an- 
alyze it by expanding it at low temperature, i.e. high 
fj. This expansion can be performed using the Euler- 
MacLaurin formula or, equivalently, a spectral represen- 
tation of the Green function G c (q,oj n ) (see Appendix 1X1 
for details). It yields 



fih (see Appendix lAlfor more details), the presence of the 
peculiar term C in (|61|) generates however odd powers 
of fih at higher orders of this semi-classical expansion. 
Indeed considering the definition of B <|23|) . one has to 
handle with care the sum over Matsubara frequencies and 
isolate the mode w„ = in the following way 



W"(Q) 

(f3h) 



y"Xi(a; m )(/Ci(a; m ) - Ki(w n + w n )) (59) 



f+OO I 

4V'"(0) / ^/Ci(w)(/Ci(w) - /Ci(a; + w„)) 



V"'(0) 2tt 
F"(0) 3 



A (g,0) + O 




with the definition 



^ofe^) = d u A Q (q, u) 

A (q,Lj) = lmG c (q,u n ^ -zw + 0+) (60) 

and where we have used the self-consistent equation for 
Io(u>n) J41|l . Although the first term in (|59|l . correspond- 
ing to the dominant one in the limit (fih) — » oo, behaves 
like w 2 and leads to a linear term, oc |cj n | in Ii(uj n ) (no- 
tice that Iq(uj +lu u ) is well defined and can be explicitly 
computed H41B 1. the second term in l|59l) is linear in u> n , 
and produces a new term C in Ii(ui n ) (we remind that 
h(uj n = 0) = by definition (gSJl ): 



/i(w n ) = (7(l-<y«, )+li(wn) 



, _ ,'ry 7T V"(o) 



(61) 



where Ji (w ra ) is a well defined function such that I\ (uj n ) oc 
for small ui n . Notice also that this term C arises only 
at finite temperature. Finally, combining Eq. (|59|l and 



> / ° g^H 



(62) 



it follows that the terms of order (T/h) 2 in (|57(l . for uj n ^ 
0, exactly cancel leading to 



Ei + h{uj n ) = 



2 V" (0) 



i + 4y"(o)/c 2 K) 
4k) 



+o 




+ 2(ld(w)-lCi(u + (j n ))j 
(63) 



Note that, although at this order Si and Ii(ui n ) admits a 
low temperature expansion involving only even powers of 



h 2 - 

B = —C 



1 



0H J q (c 9 2 + S + S 1 )(c g 2 + I]o + Ei + C) 

1 



^7(64) 



h v 

+ Jh ^ J q C q 2 + E „ + ft ( El + c) + 7 ( Wn ) + fi/i(w n ) 

Under this form, the sum over u> n can be performed 
safely, and will itself admit a low T expansion in even 
powers of fih. However, the first term in l)64|l will gen- 
erate a contribution to S 2 (and consequently also to 
-^2(^71)) proportional to (T/h) 3 , thus generating odd pow- 
ers of T. 

As we will discuss in the next subsection, (|63[) ensures 
that the temperature dependence given by the lowest or- 
der term is indeed correct. The property l|63|) can be 
shown to hold in fact to all orders^ in H and to apply also 
to other models of glasses. 



B. Specific heat 

We first focus on the expansion of (H) (|40|l to or- 
der h 2 and compute i^- The first term in 14U|) is sim- 
ply expanded by substituting S and I(ui n ) by their own 
expansion. We will leave it under this form (i.e. as 
in (|46|) by substituting So — * So + KS\ and similarly 
7o(w n ) - * Io(u n ) + fi/i(w n )), which can then be easily 
expanded at low temperature. At this order 0(h 2 ) how- 
ever, the second line of (|40|) leads to a non trivial contri- 
bution which, after some manipulations, can be written 
as 



(H) - H 



(65) 



E 



eg 2 + S + fcEi + h{un) + frh(u n ) 
fih ^ J q cq 2 + Mu> 2 + S + ftSi + I (uj n ) + hh(w„) 

V"'(0) 1 



-h 



V "(0) 2((3h) 2 ^ 



J~) Ki (ui n ) Jq (u) n )K\ (uj n ) (66) 



2h 2 i 



where -f/o, which is /37i independent is given in (|45|l . From 
this expression of (if) to order 0(h 2 ) 165(1 . we now want 
to compute the low temperature expansion to obtain the 
specific heat C V (T) to order O(h). The analysis of the 
first term in (|65() can be done similarly to ifi (|47|l . And 
as we showed previously, the finite temperature correc- 
tions to S x + I\(uo n ) behave as T 4 (|63fl (note that the 
peculiarity of the mode oj n — disappears in the first 
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term of (|65|) ~1 . one obtains that 



E 




eg 2 + S + hT,i + I (uj n ) + hli(uj n ) 
(3h ^ J q cq 2 + Mul + S + + lok) + Wi(u n ) 



(67) 



Given the complexity of the calculations, we did not in- 
tend to compute the amplitude of this term. We now 
turn to the analysis of the low temperature behavior of 
the last two terms in (I65|) . As the specific heat vanishes 
at zero temperature, the first non vanishing finite tem- 
perature correction in (|65|l is a priori of order 1 / ((3K) 2 . As 
shown in Appendix [51 it turns out that this term exactly 
cancels between these last two terms in (|65f) , leading also 
to 

-h 2 \ ^ y^/Ci(^n)y^/ (cj„)/Ci(^n) 



21r 



V"(0) 2(Ph) 2 



oc.hr - 



(68) 



Therefore, combining (|67|l and l|68[) together with (1651) . 
one obtains that 




CD ,,|) 3 + (m 5 ,^ (70) 



which shows explicitly the cancellation of the term linear 
in T in C V (T) up to order 0(h 2 ). The term fSJ gives a 
correction to the amplitude of the T 3 term in the specific 
heat gnj- 



VI. CONCLUSION 

In this paper we have studied the specific heat of a 
disordered elastic system. Using a variational approach, 
we have shown that the leading temperature dependence 
of the specific heat is C v cx T 3 at low temperatures for 
dimensions d > 2. We thus find results at variance with 
the commonly believed linear temperature dependence 
stemming from e.g. a two level system. We have com- 
puted the prefactor of the T 3 law up to the second order 
in a semi-classical expansion in K. It exhibits a depen- 
dence in the Larkin pinning length due to disorder. We 



have showed up to order h 2 that linear temperature de- 
pendence cancels in the specific heat. This property is 
quite general and is proved to all orders using a different 
approach in Refi^ where similar consequences are shown 
to hold for quantum spin glass models. 

The T 3 dependence of the specific heat also holds for 
non periodic manifolds, such as polymers or interfaces, 
whenever they can be solved by a continuous replica 
symmetry broken ansatz (or its limiting case of one step 
marginal RSB). Indeed the results we obtained in this pa- 
per do not rely on the periodicity of the elastic system. 
In the case of manifolds, the dependence of the specific 
heat we obtained becomes exact in the limit of an infinite 
number of components. 

For periodic systems, the variational method is only 
an approximation, and it would thus be interesting to 
check whether the powerlaw T 3 predicted here, has cor- 
rections away from mean-field. Such a check would need 
numerical investigations of this problem. 

For d < 2 we showed that within the framework of the 
variational method it is incorrect to use the marginal- 
ity condition criterion to compute the thermodynamic 
specific heat. The thermodynamic saddle point leads to 
an exponentially small specific heat for d < 2. As we 
discussed in this paper, this is most likely an artefact 
of the mean-field approximation. Corrections to mean- 
field should transform this gap into a pseudo-gap, with 



powerlaw density of states p(u>) 



CX U} 



leading to a spe- 



cific heat T a . The exponent a is likely to be larger than 
three in d = 1. The presence of a gap in the mean-field 
solution results from the best attempt of the variational 
method to mimic this high power. Such behavior in the 
density of states is indeed compatible with studies in one 
dimensiori 2 ^ 2 ^* 2 ^, suggesting p(ui) cx co 4 . 

An important question is thus how this density of 
states evolves when going to higher dimensions. The cal- 
culations performed in the present paper suggest that for 
d > 2 the density of states becomes p(ui) cx uj 2 . Note that 
this is at variance with the recent proposal 2 ^ that it re- 
mains p(u>) cx uj 4 as in one dimension. If this is indeed 
correct it would raise the puzzling question to under- 
stand which mechanism can lower the density of modes 
compared to mean-field. 
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APPENDIX A: EVALUATION OF MATSUBARA 
SUMS AT LOW TEMPERATURE. 



In this section, we show in detail how to extract the low 
T behavior of the sum over Matsubara frequencies in the 
equation for Ei + Ii(u> n ) (|57|l . which turns out to be cru- 
cial to understand the full structure of the saddle point 
solution. Although the standard way of studying such 
a sum is to use the spectral representation of the Green 



function, we show how to analyse it using less sophis- 
ticated method, namely the Euler Mac-Laurin formula. 
For this particular case, this method turns out to be very 
simple. Its equivalence with the conventional method is 
established in the last paragraph of this section. 



1. A first stage with Euler MacLaurin 



The aim is to evaluate the low T expansion of: 



T 



S = 27r-£«7(K|) 



(Al) 



with uj n = 2TmT/h, a Matsubara frequency, and let us 
add a cutoff |n| < N, e.g. N = hn/(2nT). We use the 
Euler-MacLaurin formula for any smooth enough func- 
tion f{x): 



N-l 



/(0) + 2 £ /(») = 2 / dkf(k) - f(N) (A2) 
j Jo 

+l(f(N) - /'(0)) - ^(f"(N) - /"'(0)) + .. 



Applying O to f(k) = 2ir(T/h)g{2kirT/h) gives: 
S-S c = 



2 L dxf{x) --A—) m 

S c = -2^g(n) + 1 {^pj g'(Sl) 



360 V h J 



(A3) 
ff"(0) 



and if the above integral converges we can take f7 — > oo 
and the boundary term S c vanishes. 



2. More sophisticated EML formula 

We now want to evaluate a more complicated sum such 
as in i|57|) : 

S = 2%—^g 1 (\2TT(n-m)T/h\)g 2 {\2TTnT/h\lA4) 

n 

2n T £ Si(|27mT/ft|)fla(|27r(n + m)T/H\) 



with y = 2irmT /h, and where we leave aside the cutoff 
from the beginning, i.e. n goes from — oo to oo. We 
choose to > but similar calculations hold also for m < 
0. We use the decomposition: 



S = Si + s 2 + s 3 



(A5) 

+2irj gi (\27rmT/h\)g 2 (0) + 2irj gi (0)g 2 {\2TrmT/h\) 
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with 



Sl=27: T gi(\2n(n-m)T/h\)g 2 (\2nnT/h\) 



S 2 = 2nj J2 9i(\Mn - m)T/h\)g 2 (\2imT/h\) 



n— — oo 
m — 1 



S z = 27T T E 9i(\2n(n-m)T/h\)g 2 (\2TrnT/h\) 



n=l 
00 



n— m+1 



Applying the standard Euler-MacLaurin formula <A2() to 
Si one obtains up to terms of order 0((T//i 4 )) 



T 



Sl = 2n -fr H -9i( 27r ( m + n)T/h)g 2 {2irnT/h) 

n=l 

00 \ T 

dxgx{x + y)g 2 {x) - -^27r—g 1 (y)g 2 (Q) 

+(2irj) 2 -^(-g[(y)g 2 (0) - gi{y)g'M) (A6) 



similarly to S 2 



S 2 = 2tt% V 5i(2tt(to - n)T/h)g 2 (2irnT/h) 



rv it 

dxgi{y - x)g 2 (x) - -2tt— gi(y)g 2 (0) 
io z n 

-l2n^ gi (0)g 2 {y) 

-(2^^) 2 ^(5i(y)ff 2 (0) - ff!(y)52(0) - g'x(0)g 2 (y) 



-gMg'M) 



(A7) 



and 5s 



T 



S3 = 2tt- ^ 5i(27r(n - m)T/K)g 2 {2imT/h) 

n— m+1 

00 IX" 

c?X5i(x - y)g 2 (a;) - ^2Tr-g 1 (0)g 2 (y) 

+(2nj) 2 ±(-g[(0)g 2 (y) - g 1 (0)g' 2 (y)) (A8) 
Collecting the terms in 1A6I lATI !a"8)) yields finally 



S= / dxgida; - y|)p 2 (a;) + / dxgx{x + y)g 2 (x) 
Jo Jo 

-2(2^|) 2 l(. 9l (y)^(0) + <?i(0) 52 (y)) + 0{{T/hf) 

(A9) 



3. Euler Mc-Laurin vs spectral representation 

We want to compute the high j3 expansion of the fol- 
lowing Matsubara sum which enters the equation l|57|) : 

§(u n ) = E ^1 (w m )/Ci (w„ +w m ) (A10) 



The Euler-MacLaurin formula established previously 
(|A9I) allows us to compute the first term of this expan- 
sion, namely of order (T/H) 2 : 



S(uJ n ) = S(ui n )\/3h=c 

1 2tt 

+ 



(All) 



m)2 3 -^/o(x)U^ 2 ( S o)/C 1 K) + 0( w ) 
T\ 2 2ir 



S(ui n )\f3h 

— OO ~T 



h 3 



4,(0,g)£i(w n ) 



+0((T/S) 4 ) 

where ^4o(w, g) is given in the text l|60() . We want to show 
that this term of order (T/H) 2 in IjAlljl can be obtained 
in a more standard way using a spectral representation 
to compute the sum over Matsubara frequencies. Indeed, 
we can write this sum (|A10|) as 

S{u n ) = I { — / — A {q,ui)A Q (q , ,u 2 ) 

J q : q' J — OO " J —oq " 

' (A12) 



/3ft ^— ' zcj m — ui zw m + iu> n — u 2 

TCI 

with Ao(q, uj) given in (|6(J|) . The sum over the Matsubara 
frequencies is then straightforwardly computed, and it 
gives 

— / —A (q,ui)A a (q',u 2 ) 



S(u n ) = 

f q,q' J — oo 

/s(iti) - Jb{u 2 



(A13) 



U 2 — U\ — lLV n 

dui I du 2 —A (q 7 u 1 )A (q',u 2 ) 

> q,q' J — 00 J — 00 ^ 

x/b(«i) ( 1 — + 1 —— ) (A14) 

\ U 2 — U\ — lU! n U 2 — U\ + lOJ n J 

The term in (T/ft) 2 in S(u> n ) is obtained by computing 
j^TTgfp 1 where the derivative concerns the explicit de- 
pendence in /3ft, i.e. does not act on the implicit one of 
uj n . This leads to 



ai/(/?n) 2 

(/3ft) 3 



(A15) 



1 



. } , , du! / du 2 —A (q,ui)A (q',u 2 ) 

■' q,q' J — 00 J — 00 " 



1 



1 



"2 — Ui — iuj n ' u 2 - ui + iuj n J Ul (eP hui - l) 2 

The integral over u\ is well behaved due to the derivative 
of the Bose factor and we can safely rescale u\ — > (3hxi : 



dS(uj n _ 

di/(ph) 2 ~ 

p p r+00 r+00 

2 



(A16) 



q,q' J —00 



dxi / du 2 —A (q,xi/ f3)A (q' ,u 2 ) 



1 



1 



u 2 - x\/f3 — iuj„ u 2 — x\ J (3 + iuj n ) (e xi - 1) 



xi- 
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From that expression, we extract the coefficient of order 
l/{(3h) 2 : 



<91/(/3ft) 2 
1 



(3K — >OC 



4»(?,o) 



1 



+ OC 



' q j — oo 

i 



du 2 



A (q,u 2 ) 



U 2 — lLO n U2 + lUJ n 



2tt 



2*1 - l) 2 

(A17) 



where we have used the value of the integral 



dxx £ 



(e x - l) 5 



This calculation shows explicitly the equivalence 
(IA11I) . I|A17|) of the two methods to compute this low T 
expansion. The expansion of the term entering the equa- 
tion for Si + Ix(u> n ) can be written 



4V'"(0) 



08ft) 

= 4V 7 " (0) 



/Ci (u m ) (fCx (w m ) - /Ci (u n + o; n )) (A19) 
' 50 dw. 



2tt 



/Ci(w)(/Ci(w) - /Ci(w + w„)) 



(£) V'"(0)^^^( (? ,0)( t 7 1 (So)-/C 1 K)) 



+0((T/h) 4 ) 



(A20) 



And using the equation for I (uj n ) given in the text (|41ll . 
the term of order (T/h) 2 in (|A19|) can simply be written 



4V"'(0) 



(/3ft) 

= 4 V"' (0) 



^/Ci(w m )(/Ci(w m ) - /Ci(w„ + w„)) (A21) 

m 



2tt 



/Ci(w)(/Ci(w) - /Ci(w + w„)) 



as given in the text l(5§fl. 



APPENDIX B: LOW TEMPERATURE 
EXPANSION TO ORDER C»(ft 2 ) : DETAILED 
CALCULATIONS. 

In this appendix, we focus on the internal energy to 
order 0(h 2 ) H65JI and show how to extract the coefficient 
of the term cx \/(j3K) 2 in the expression 



V"(0) 2{(3h) 2 ^ 
2ft 2 1 ^ 



^1 (^n ) ^2 h (W„ )/C X (W„) 



(/3ft) 2 



(Bl) 



Let us expand the first term of ijBTJ: 
V"'(0) 1 



V"(0) 2(/?/i) 2 



Y, Ki K) J] *o(w»)£i K) (B2) 



The general structure of the high (3h expansion of this 
term is the following 



T(2) 1 

1 = X (0) + — + Of—) 



(B3) 



We are not interested in the constant as it does not con- 
tribute to the specific heat and and focus here on the 
computation of the first term X 1 - 2 ' . We use a spectral 
representation to compute the Matsubara sums in 2 : 



1 



c^ + Eo + Mw 2 +/o(w n ) 



_1 r+°° i 
— / duA (q,u)- 



/Ci(w„)/ (cj n ) 



cq 2 + E + M W 2 + J ( Wn ) 



(B4) 



(B5) 



-1 r°° 1 
= / du{I%(u)B (q, u) + %(u)Ao(q, u))- 

7T J_ OQ lLJ n - U 

where Ao(q,uj) is defined in the text l|60|l and 

Bo(q, oj) = KeG c0 (iuj m — > u> + iS) (B6) 
c9 2 -w 2 + S + /qM 



(aZ 2 -u, 2 + £ + J6(u)) 2 + /£W 2 
Using the identity 

n 

X can be written 
-V"'(0) 1 



(B7) 



(B8) 



1 = 



V"(0) 2. 2 y_ dmJM^)f B (u x ) (B9) 



x / du 2 (Io(u 2 )B Q (q,u 2 ) + Io(u 2 )A (q,u 2 ))fB(u2) 

J — oo J q 

We use this expression to compute the term of order 
l/(/3ft) 2 , together with dl/ ^ h) 2 = -^^gpij' we have 



01 



<91/(/3ft) 2 

-V'"(0) 1 (/%) 3 '•+ 00 



(BIO) 



v^W^—J^ dMl i/ o(g ' Ml)ui ( e ^-i) 2 



a f3hu 1 



+00 



x / du 2 (Ig(u 2 )B Q (q,u 2 ) + lQ(u 2 )A (q,u 2 ))f B (u 2 ) 



V"(0) 2 

x y (Io(u2)B (q, u 2 ) + ^( M2 )A (g, u 2 ))M2 _ ^ 2 
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After some manipulations, we obtain this coefficient : 
dl 



(3h=oo 



dl/(l3h) 2 

x / (7£(ui)B (g,ui) +/o(«i)A)(?.«i)) 



(Bll) 



+ / (0Jft)(O)B o (O,g) / du 2 / A)(3,«2)) 

J q J — oo «/ g 



with a given in (|A18() . We can simplify that expression 
using the equations for Iq(u) and Iq (u) obtained from 
JUJ: 



l£(u;) = 4K"(0) / A (q,uj) 
'q 



(B12) 



which finally yields 
AT 



3l/(/3fi) 2 



fjh—oo 



(B13) 



y"'(o) a 

7"(0) 2tt 



2 (a u /")(0) / / duA (q,u)B (q',u) 

q.q 1 J — oo 



The high f3h expansion of the following term 



J = \v'" (o) Y KiK)^! K)^! K + w «) 



"'I'm /-+ 00 



21/'" (0) 
3tt 3 



— Y — 



du 1 du 2 du 3 A (u 1 )A (u2)A (u 3 ) 



(B14) 



is performed using the same kind of manipulations, which 
lead to 



dj 

I' (u) = -iV"(0) ( [ 2 \ - B (q,cj)) dl/m 

\Jq -h Zjq / 



(3h—oo 



V"'(0) a 
V"(Q) 2tt 2 



MM 



(B15) 



duA (q,u)B (q' ,u) 



q,q' J — oo 



Combining (|B13|I and (|B15|) shows that the term of order 
{T/hf cancels in (|BlTl . as given m the text (jHEJ- 



